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POTENTIAL AND VISCOUS FLOW PREDICTION IN 
V/STOL PROPULSION SYSTEM INLETS 

by 

Norbert 0. Stockman 
Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 


ABSTRACT 

Highlights of the method of analysis of inlet flow are given. To 
indicate the accuracy of the method, several comparisons with experiment 
for different V/STOL inlet configurations and various operating condi- 
tions are given. Two applications to inlet design and analysis are then 
discussed. A summary of current efforts is given, and finally areas of 
possible future work are indicated. 
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INTRODUCTION 

Over the past several years a method has evolved at the Lewis Re- 
search Center for the calculation of potential and viscous flow In 
subsonic inlets at arbitrary operating conditions. The development 
of tho potential flow calculations (ref. 1) was originally motivated 
by the need to design efficient inlets for an in-house lift fan test 
program (refs. 2 & 3). The method was quite successful at this and was 
eventually extended to several other VTOL applications (ref. 4). 

When the method was applied to STOL inlet designs and analysis, 
the boundary layer became important, and the boundary layer (or viscous) 
calculation (ref. 5) was incorporated into the procedure." A status 
report on the resulting method as of late 1973 is given in reference 6. 
Since that report many additional applications were made (e.g., refs. 

7 & 8). 

The present paper will present the highlights of the method itself, 
several comparisons with experiment, some recent applications of inter- 
est, a summary of current efforts and finally a discussion of possible 
future work. 


SYMBOLS 


A, B, C 
a, b 

C. 


D 

M 

P 

P, Q 
s 

V 

w 

a 

6 * 

f 


combination coefficients, eq (1) 

semi-major and minor axis of bisuperellipse, eq (3) 
skin friction coefficient (ratio of wall shear stress 
to dynamic pressure at edge of boundary layer) 
diameter 
Mach number 
pressure 

bisuperellipse exponents, eq (3) 

surface distance 

velocity 

Inlet mass flow 

inlet incidence angle or angle of attack 
boundary layer displacement thickness 
density 


c 

cor 

h 

i 

J 

ref 


t 

th 


00 


Subscripts: 

control station 

corrected for compressibility 

highlight 

incompressible 

basic solution, j ■ 1, 2, 3 

reference value 

total or stagnation conditions 

throat 

frees tream value 
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Superscripts: 


average value 
vector quantity 


METHOD OF ANALYSIS 

The method of analysis is well documented in the literature (refs. 
1, 4, 5, 9 & 10), and only the highlights will be presented here. The 
method comprises the following several steps (shown schematically in 
fig. 1): 

1. Geometry representation (SCIRCL in fig. 1) 

2. Incompressible potential flow basic solutions (EOD) 

3. Combined solutions with compressibility correction (COMBYN) 

4. Boundary layer calculations (VISCUS) 

5. Iterative loop 

The four computer programs are available from COSMIC, Computer 
Center, Information Services, 112 Barrow Hall, University of Georgia, 
Athens, Georgia, 30602. Programs SCIRCL, EOD and COMBYN are one unit 
with number LEW-12152; program VISCUS is number LEW- 12178. 

Geometry 

The inlet is assumed to be axisymmetric and is represented by its 
meridional profile. This profile is broken into segments at convenient 
tangent points as shown in figure 2. Each segment may be defined by an 
analytic expression or a set of points. The inlet duct walls and the 
outer surface (nacelle or belimouth) must be extended far downstream 
(fig. 2) to facilitate obtaining accurate potential flow solutions in 
the inlet region of interest. The geometry program (SCIRCL, fig. 1) 
prepares coordinate-point input for the potential flow program and also 
prints out information such as curvature, wall angles, area distribu- 
tion, etc., which is useful in preliminary screening of proposed inlet 
shapes. 


Incompressible Potential Flow 

The Douglas axisymmetric incompressible potential flow computer 
program (EOD, ref. 9) is then used to obtain three, independent basic 
solutions. These three basic solutions V,, j = 1, 2, 3 are then com- 
bined (program COMBYN, fig. 1) into a solution of interest having arbi- 
trary flow conditions of V ffi , a, and mass flow W (fig. 3). The 
combination equation is: 
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V - av l + bv 2 + riT (l) 

where A, B, and C are obtained from the flow conditions. Thus, once 
the basic flow solutions are obtained for a specified geometry, any 
solution of interest for that geometry can be obtained without repeat- 
ing the more time-consuming potential flow calculations. 

Compressibility 

by equation (1) is incompressible apd ia cor- 
by the Lieblein-Stockroan compressibility 


_ V / V 

v i ( V oi 1 (2) 

right hand side are obtained from the incom- 
the input flow conditions. From the velo- 
(Mach number, pressure ratio, streamlines, 


Viscous Flow 

The surface Mach number distributions obtained from program 
COMBYN are used as input to the Herring-Mel lor boundary layer cal- • 
culation (program VISCUS, fig. 1). Reference 5 contains a complete 
documentation of program VISCUS and references to the original sources. 
Program VISCUS calculates boundary layer profiles, displacement thick- 
ness 5 *, skin friction coefficient Cf, etc., at each station, and alao 
predicts transition from laminar to turbulent flow and separation 
(whether laminar or turbulent). 

Iterative Loop 

If the boundary layer is significant in the inlet region of in- 
terest it may be desireable to add the displacement thickness 5 * to 
the original inlet profile and repeat the entire' solution procedure, 
thus obtaining a new Mach number distribution, new 6*, etc.. Thin 
process may be iterated to satisfactory convergence. Usually one 
iteration is sufficient. In parametric studies or preliminary design 
screening often no iteration is needed. 

COMPARISON WITH EXPERIMENT 

To indicate the accuracy of the various aspects of the prediction 
method several comparisons with experimental data will be given. 

These will range over incompressible flow, compressible flow with 


The velocity obtained 
rected for compressibility 
correction (ref. 11). 


V 

cor 

where all the terms on the 
pressible flow solution or 
city other flow properties 
etc.) are obtained. 
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insignificant boundary layer and viscous flow with various conditions 
of the boundary layer. 


Incompressible 

The surface pressure distributions in a chordwise cut of a fan- 
in-wing inlet ore shown in figure 4. Three surfaces are Bhown on the 
plot; the forward surface of the bellmouth, the centerbody and the aft 
surface of the bcllmouth. The agreement between theory and experiment 
is quite good everywhere on the inlet. This case is included to illus- 
trate the adequacy of the method of geometry representation even when 
the calculational model (see'fig. 2(b)) differs rather significantly 
from the real geometry. 


Compressible 

The next case illustrates the applicability of the method when 
the flow is compressible. Figure 5 shows the theoretical and experi- 
mental surface pressure distributions on a fan-in-pod inlet (described 
in ref. 3). Both incompressible and incompressible-corrected-for- 
compressibility theoretical curves are given. The experimental static 
pressures agree quite well with the theory corrected for compressi- 
bility along the entire surface of the inlet. 


Viscous 

A translating-centerbody sonic inlet in the retracted configuration 
(fig. 6(a)) is used for comparing theoretical boundary layer results 
with experimental. Various conditions of the boundary layer are ob- 
tained by varying the inlet incidence angle a. In all cases, the one- 
dimensional throat Mach number iB 0.50 and the free stream Mach number 
is 0.13 (these results are all taken from ref. 8 ). 

Boundary layer attached . - Surface Mach number distributions for 
zero incidence angle are shown in figure 6(b). Two theoretical curves 
are shown to illustrate the effect of the boundary layer on the surface 
Mach number. The solid curve is the potential flow alone (no cor- 
rection); the broken curve is the potential flow obtained with 8 * added 
to the inlet profile and is in excellent agreement with the theory. 

No separation is indicated either by theory or by experiment. 

Diffuser separation . - Results are shown in figure 6(c) for an in- 
cidence angle of 40°. Here the theory, with or without the 6* correc- 
tion predicts separation at about the same point as the experimental 
data indicate the start of a "pressure plateau" commonly associated 
with observed separation. 


Lip separation , - On figure 6(d) results are given for an inlet 
incidence angle of 50°. The theoretical Mach number distribution i* 
what would be obtained if the inlet did not separate. The experimental 
data clearly indicate separation on the inlet lip. The theory also 
predicts lip separation and although possibly not at the exact point 
as experiment, it is probably adequate as a guide in inlet design. 

. In cases where the theory predicts separation the calculations 
stop and there is no 6* beyond the point of separation. In such cases 
the distribution of 6* into the separated region is obtained by extra- 
polation using an unseparated case (lower a) as a guide ^see ref. 8). 

Additional comparison for this and other configurations can be 
found in reference 8. On the basis of these results, the theory seems 
adequate in predicting the boundary layer behavior for this type of 
inlet configuration. 


RECENT APPLICATIONS 

Two recent applications of tha theoretical method to inlet design 
and analysis will be discussed. One involves only potential flow and 
the other viscous flow. 


Potential Flow 

This example is taken from a theoretical screening study of in- 
lets for the Quiet Clean Short-Haul Experimental Engine (QCSEE) Project 
(ref. 7). This study was undertaken to ensure that a reasonable ,‘ange 
of contraction ratios was chosen for experimental evaluation. Figure 
7 shows the effect of the internal lip contraction ratio, (D h /D t ), on 
the surface Mach number distribution at QCSEE takeoff conditions (M^" 
0.79, Mo,, = 0.12, a - 50°). The figure shows that the higher the con- 
traction ratio, the better the Mach number distribution, i.e., the 
lower the peak and the less severe the unfavorable gradient. (However, 
the takeoff performance must be compromised with the cruise perfor- 
mance, and the highest contraction ratio lip may not be beBt overall.) 


Viscous Flew 

The sonic inlet of figure 6(a) will be used to illustrate an 
application of the viscous flow calculation to inlet design. Tha 
distribution of theoretical skin friction coefficient as a function 
of surface distance from stagnation point to diffuser exit is shown 
in figure 8 (taken from ref. 8). The flow conditions are - 0.50 
and Moo“ 0.13 at several values of incidence angle. 
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The criterion for separation prediction is that the Bkin friction 
go to zero, therefore a local minimum with a low value of C£ is a point 
of likely separation as the operating conditions become more severe. 

The zero incidence curve indicates two regions of possible boundary 
layer separation (at the two minimums), one in the diffuser (where tha 
flow is turbulent) and the other on the internal lip (where the flew 
is initially laminar). 

At an incidence angle of 20°, separation is indicated in the dif- 
fuser. As incidence angle is increased further up to 40° the separ- 
ation point moves upstream somewhat but remains in the diffuser at a 
position close to the diffuser wall inflection point (S ~ 0,6). At 
50° incidence, however, th« separation point jumps to the lip. 

Experimental results for this inlet indicate that diffuser separ- 
ation in the region of S *» 0.6 leads to only small losses whereas lip 
separation leads to intolerable losses. Thus, it appears that a dis- 
tribution of C^, like that of figure 8, is undesireable and if possible 
the inlet design should be modified to raise or eliminate the minimum 
on the lip. 


CURRENT WORK 

This section will briefly discuss several areas of ongoing work 
and work planned for the immediate future. These areas fall under two 
general headings: investigation of the method itself; and application! 

of the method. 


Investigation of Method 

Potential flow . - Improved methods of solution of the incomprea- 
sible potential flow were recently obtained on contract from Dougla! 
Aircraft Company (ref. 12). These improvements were incorporated into 
the calculation procedures and are currently being evaluated both by 
internal program checks and by comparison with experiment. Preliminary 
results indicate increased accuracy and reduced computer times. 

Viscous flow . - The flow in an axisymmetric inlet at angle of 
attack is three-dimensional, that is, in general, there is a circum- 
ferential component. The boundary layer calculations being two- 
dimensional neglect this circumferential component. Along the windward 
meridian there is no potential flow circumferential component and the 
method should be at its best here. The data presented herein seem to 
indicate that the method is adequate here. At other circumferential 
locations, however, the method may be inadequate since the boundary 
layer calculations proceed along a meridian and neglect the circumfer- 
ential component of flow and the consequent spiral streamline. 
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To evaluate the adequacy of the method at various meridians, an 
investigation is underway in which experimentally-determined boundary- 
layer profiles obtained at several circumferential positions will be 
compared with calculated profiles. Hopefully this study will establish 
conditions under which the method is adequate and make recommendation* 
for changes in the method where it is not adequate. 

Applications 

Scale effects . - Host inlet testing, particularly that involving 
setting separation bounds, is done on relatively small scale model*. 

It is generally not feasible to adequately investigate scale effect* 
experimentally, certainly not up to full scale. However, scale effect* 
can be investigated to a limited degree by the existing program. Th* 
effects of Reynolds number and wall curvature are included in the pro- 
gram, and thus geometrically similar inlets at different scale will 
have different calculated boundary layer development. At present th* 
program does not include the effect of scale on shock-boundary layer 
interaction or on more complex transition mechanisms (such as laminar 
separation with turbulent reattachment). An investigation is now 
underway to study the effects of scale within the present limitation* 
and to recommend improvements in the modeling of the transition pro- 
cess. It is hoped that eventually the boundary layer behavior at any 
scale may be confidently predicted. 

Unsymmetrical inlets . - Some inlets of current interest are not 
axisymmetric, e.g., the scoop inlet (shown in figure 9) designed for 
reduced downward directed noise and low foreign object damage. A 
study is getting underway to investigate means of approximately ana- 
lyzing the flow in such inlets using the axisymmetric programs. If 
the skew angle is relatively small an equivalent angle of attack 
approach will be used, i.e., the inlet will be analyzed as if it 
were an axisymmetric inlet operating at sane equivalent artificial 
angle of attack related to the actual angle of attack and the inlet 
skew angle. 

Tilt-nacelle lip design . - Some V/STOL propulsion system design* 
may involve a tilting nacelle. This means the inlets must operate 
efficiently under the conflicting requirements of vertical flight 
(Vm = 0), crossflow (with conditions at transition possibly as severe 
as = 120 knots and a «= 100°) and at high speed cruise. The fea- 

sibility of using a fixed geometry lip on this type of inlet will be 
studied by theoretically investigating several possible lip shapes. 

Bisuperellipse diffuser walls . - Most of the diffuser wallB used 
in experimental and analytical studies at Lewis have been cubics. 

Once the diffuser length and the fan and throat diameter have been 
chosen a cubic is determined. However, with a bisuperellipse 
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(x/a) P + ,(y/b) Q - 1.0 (3) 

only a and b are determined and P and Q are still free. This means 
that the location of the inflection point and its slope (i.c., the 
maximum wall angle) can be specified. (When either P or Q is lest 
than 1.0 and the other greater than 1.0 the curve has an inflection.) 
It is proposed to do a parametric study of the effect of the location 
of the inflection point and the maximum wall angle with the hope of 
obtaining shorter and more efficient diffusers than those in current 
use. 


FUTURE POSSIBILITIES 

A few areas of possible future study will be briefly mentioned. 
These studies would be made only if a real need is indicated and if 
sufficient time and manpower are available. 

Modification of Potential Flow Calculations 

It seems feasible to modify the calculations in the COMBYN pro- 
gram to account in an approximate way for local total pressure losses 
along the walls and possibly also the effect of shocks. 

Modification of Boundary Layer Calculations 

As was stated previously the boundary layer calculation precedes 
along a meridian and neglects any circumferential velocity component. 
However, this circumferential velocity component is available in the 
potential flow results and could be used to construct a spiraling 
three-dimensional streamline along the inlet surface. The velocity 
distribution along this streamline could be input to the boundary 
layer program. There the boundary layer calculation could proceda 
along a more realistic streamline and the results should better approx- 
imate the real flow, at least to the degree that the current method 
approximates the flow along the windward meridian. 

Another shortcoming of the boundary layer calculation is the 
neglect of shock interactions. Since many cases of current interest 
contain regions of local supersonic flow it may be necessary to account 
for possible shock-boundary- layer interaction in the boundary layer 
calculations • 

Finally, the program is not able to predict separation bubbles 
that appear to be present in the experimental data. However it might 
be feasible to include the predictions of a separation bubble in an 
improved transition model. 
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Three-Dimensional Analysis 

It seems likely that eventually a more accurate analysis of unsym- 
metrical inlets will be required than that obtainable from the axl- 
symmotric programs. In this eventuality, the Douglas three-dimensional 
potential flow program could be modified to handle the desired arbi- 
trary input flow conditions (V ra , a, $) and to include compressibility 
and boundary layer calculations comparable to those now included in tht 
axisymmetrlc calculations. 


CONCLUDING REMARKS 

• 

A method for the prediction of potential and viscous flow in axi- 
symmetrlc propulsion system inlets has been briefly described. The 
method has already proven to be a very useful and powerful tool for 
both analysis and design purposes. The various elements of the method 
are constantly being updated, thus improving both its versatility and 
accuracy. It is probably the best tool available at present for 
analyzing the compressible viscous flow in axisymmetrlc inlets at 
arbitrary angle of attack and will continue to be until exact three- 
dimensional compressible inlet flow computer programs become practical. 
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Figure 2 . - Representative inlet gwmetries lor potential Mow solutions 
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figure 7. - filed a I internal lip contraction ratio on surface Mach 
number distribution. 



Figure 8. - Theoretical skin friction coefficients 
(without 6' correction for centerbody retracted. 
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